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Abstract 

We introduce operator definitions for virtual photon fragmentation functions, 
which are needed for reliable calculations of Drell-Yan transverse momentum 
(Qt) distributions when Qt is much larger than the invariant mass Q. We 
derive the evolution equations for these fragmentation functions. We calculate 
the leading order evolution kernels for partons to fragment into a unpolarized 
as well as a polarized virtual photon. We find that fragmentation functions 
to a longitudinally polarized virtual photon are most important at small z, 
and the fragmentation functions to a transversely polarized virtual photon 
dominate the large z region. We discuss the implications of this finding to 
the J/ip mesons' polarization at large transverse momentum. 

PACS Numbers: 12.38.Bx, 12.38.cy, 13.85.Qk, 14.70.Bh 
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I. INTRODUCTION 



Gluon distribution plays a central role in calculating many important signatures at 
hadron colliders because of the dominance of gluon initiated subprocesses. A precise knowl- 
edge of the gluon distribution is absolutely vital for reliable predictions for signal as well 
as background cross sections |JJ. A great effort has been devoted to find good physical 
observables for extracting information on gluon distribution 

For many years, prompt photon production has been thought as a clean signal for infor- 
mation on gluon distribution because its cross section at lowest order is dominated by the 
"Compton" subprocess: g + q — > 7 + q, and this dominance is preserved at higher orders 
However, the theoretical and experimental complications have limited our ability to 
extract clean information on gluon distribution from the direct photon data j5|. At col- 
lider energies, prompt photons are observed and their cross sections are measured only if 
the photons are relatively isolated in phase space. Isolation is required to reduce various 
hadronic backgrounds. But at the same time, the cross section is no longer totally inclusive 
and theoretical predictions become sensitive to the isolation parameters f|,[|. In addi- 
tion, phenomenological fragmentation functions are needed for including photons emerged 
from the long-distance fragmentation of quarks and gluons that are themselves produced 
at short-distance 0,0]. Our knowledge on phenomenological fragmentation functions and 
the theoretical uncertainties associated with the isolation prevent fully quantitative deter- 
minations of gluon distribution from the collider data on isolated photons. Although data 



at fixed target energies provide good information on gluon distribution at large x | 10[] , the 
controversy about how much /^-smearing is required to understand the data introduces the 
significant uncertainties to the gluon distributions |TT| , p^| . As a result, all direct photon 



data were excluded from recent CTEQ global analyses of parton distributions H . 

Recently, Berger, Gordon, and Klasen (BGK) show that Drell-Yan transverse momentum 
{Qt) distributions in hadronic collisions are dominated by partonic subprocesses initiated 
by incident gluons if Qt > Q/2, where Q is the invariant mass of the produced lepton pairs 
|13|| . BGK argue that Drell-Yan Qt distribution is an advantageous source of constraints on 



the gluon distribution, free from the experimental and theoretical complications of photon 
isolation that beset studies of prompt photon production. 

Other than the difference between a virtual and a real photon, Drell-Yan and prompt 
photon production share the same partonic subprocesses. The virtual photon produced in 
Drell-Yan process subsequently decays into a pair of leptons. Since we are mainly interested 
in the cross section at high Qt and low Q, we ignore the Z channel contributions in this 
paper. If we integrate over angular dependence of the lepton pairs, the Drell-Yan massive 
lepton-pair production between hadron A and B can be expressed in terms of an inclusive 
production of a virtual photon [|13[ 



d(JAB^l+l-{Q) _ ( OLem \ d(?AB^*{Q) 

dQ 2 dQ 2 T dy ~ \3irQ 2 j dQ\dy 



(1) 



Because of the advantage of measuring the leptons, Drell-Yan massive lepton-pair production 
as well as the inclusive virtual photon production defined in Eq. ([I]) are entirely inclusive. 
The usual factorization theorems in Quantum Chromodynamics (QCD) should apply |l4l , |T5 
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d(T ^Q2 J d y Q) = J2 J dX!(j) a/A (Xi, H) J dx 2 <k/B(X2, fl) da ^ 1 ^y ) X 2, Q, Qt, V\ fl) (2) 

where J2ab run over a U parton flavors, the <fi a /A and <pb/B are parton distributions, and 
/i represents both renormalization and factorization scale. In Eq. (|2[), da a b^*(Q)/ ' dQ\dy 
are short-distance partonic hard parts and perturbatively calculable to all orders in a s (n). 
Similar to prompt photon production, the lowest order "Compton" subprocess to a virtual 
photon: g + q — > -y* + q dominates the Qt distributions at large Qt as long as the collision en- 
ergy is high enough to overcome the phase space penalty caused by the virtual photon mass. 
Therefore, Drell-Yan Qt distribution at large Qt is an advantageous source of information 
on the gluon distribution |T3|| . 

When Qt 3> Q, the Qt distributions calculated order- by-order in a s in the conventional 
fixed-order QCD perturbation theory receive a large logarithm, ln(Q T /Q 2 ), at every power 
of a s beyond the leading order. Therefore, at sufficiently large Qt and \/S, the convergence 
of the conventional perturbative expansion in powers of a s is impaired, and the logarithms 
must be resummed. 

In order to resum the large logarithm, we introduce a concept of virtual photon fragmen- 
tation functions Df_+ 7 *(z, fip] Q) for a parton of flavor / to fragment into a virtual photon 
of invariant mass Q. Normally, a virtual particle state is not physical, and therefore, a frag- 
mentation function to such a state may be gauge dependent and ill-defined. However, if such 
a virtual state immediately decays into a completely measured physical state, we believe that 
a fragmentation function to such a virtual state is effectively physical. The fragmentation 
function is experimentally measurable if the decay to the physical state is calculable. 

Unlike the real photon fragmentation functions [[|, the virtual photon fragmentation 
functions are fully perturbative if Q 3> Aqcd- In terms of the virtual photon fragmentation 
functions, the conventional perturbative expansion for d<7 a &_,. 7 *(Q)/ dQ T dy in Eq. (^j can 
be reorganized according to a new factorization formula such that the large logarithms are 
resummed to all orders in a s . The detailed derivation of the new factorization formula for the 
virtual photon production at Qt > Q and the discussions on the all order resummations will 
be published elsewhere. In this paper, we concentrate on the process independent physics 
associated with the virtual photon fragmentation functions. 

In the next section, we derive the cut-vertex and corresponding operator definitions for 
virtual photon fragmentation functions. From the definitions, we derive the evolution equa- 
tions (or renormalization group equations) for these fragmentation functions. We calculate 
the leading order evolution kernels for the evolution equations. By solving the evolution 
equations with the calculated evolution kernels, we derive the virtual photon fragmenta- 
tion functions. We show that the virtual photon fragmentation functions are in principle 
perturbatively calculable to all orders in a s if Q ^> Aqcd- 

Because of our ability of measuring the leptons, we can probe the polarization of the 
virtual photon in Drell-Yan massive lepton-pair production. Therefore, it is also physically 
meaningful to define fragmentation functions to a virtual photon with a specific polarization. 
In Sec. PTT| , we derive the evolution equations and the leading order evolution kernels to 
polarized virtual photon fragmentation functions. 

In Sec. [TV], we present our numerical results for virtual photon fragmentation functions. 
By showing the fragmentation functions at different scales, we demonstrate the evolution 
properties of the fragmentation functions for partons to a unpolarized as well as a polar- 
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ized virtual photon. We find that fragmentation functions to a longitudinally polarized 
virtual photon are most important at small z, and the fragmentation functions to a trans- 
versely polarized virtual photon dominate the large z region. When Qt is large while \/~S 
is fixed, fragmentation functions at large z are more relevant for calculating the cross sec- 
tions. Therefore, we conclude that the virtual photons produced in a unpolarized Drell-Yan 
massive lepton-pair production are more likely to be transversely polarized at high Qt- 

Recent data on J/ip polarization measured by CDF collaboration at Fermilab Tevatron 
seem to be inconsistent with the predictions from various models of J/i[) production ||16|| . 
The virtual photon production (extracted from Drell-Yan massive lepton-pair production) 
at large Qt and small Q 2 has a lot in common with the J/ip production at high Qt- They 
both have two large physical scales: Qt and Q 2 , which is equal to M 2 ^ in the case of J/ip 

production; and Q T is much larger than Q 2 . If the collision energy v^S" is large enough and 
the logarithm \n{Q T /Q 2 ) is so important that the resummed fragmentation contributions 
dominate the production cross sections, the virtual photon and J/i/j production will share the 
same partonic subprocesses. Only difference between the virtual photon and J/i/j production 
at high Qt is the fragmentation functions. The virtual photon fragmentation functions are 
completely perturbative, while the parton to J/ip fragmentation functions involve final-state 
nonperturbative interactions. Understanding the difference in such final-state interactions 
is very important for reliable predictions of J/i/j production. We propose to measure the 
virtual photon polarization in Drell-Yan massive lepton-pair production at large Qt and low 
Q 2 . Because the virtual photon polarizations in Drell-Yan massive lepton-pair production 
are completely calculable and independent of the final-state nonperturbative effect, the 
measurements of the virtual photon polarizations at high Qt provide not only a good test of 
QCD perturbation theory, but also a reference process to test the models of J/ip formation. 



II. UNPOLARIZED VIRTUAL PHOTON FRAGMENTATION FUNCTIONS 

Like other fragmentation functions |T7|], a virtual photon fragmentation function, 
Df^*(z, fiF',Q) is defined as a probability density to find a virtual photon of invariant 
mass Q and momentum fraction z from a parent parton of flavor /. In this section, we 
derive the operator definitions for virtual photon fragmentation functions, and correspond- 
ing evolution equations. We calculate the leading order evolution kernels for the evolution 
equations and derive the virtual photon fragmentation functions by solving the evolution 
equations numerically. 

In order to simply our derivations, we choose a frame in which the virtual photon is 
moving very fast along the z-axis, 

Q" = (Q+Q- T ), and Q~ = ^ (3) 

with Q + ^> Q~ . We also introduce two useful vectors 

n" = (l,0,0r) and r^ = (0,l,0 T ) (4) 
with n 2 = n 2 = and n^n^ = 1. For any four- vector p, we have p M n M = p~ and p^n^ = p + . 
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Consider a generic quark to virtual photon fragmentation process, as shown in Fig. |T](a). 
The top part, labeled by T, corresponds to the fragmentation from a quark of momentum 
A; to a virtual photon of invariant mass Q; and the bottom part, labeled by B, represents a 
short-distance hard collision at a scale fi ^> Q. By carrying out collinear expansion of the 
quark momentum k in the B at k^ = (Q + /z)n 11 and separation of the trace between the top 
and the bottom we factorize the generic fragmentation process as follows, 



d 4 k 

W) 



4 Tr [B g (f,,k)T q ^(k,Q)) 



— Tr 

z 2 



B q (^k = ^)\ 1 -l^- 



d A k 



Tr 



z 



4^ 



z 2 S(z- 



z 

Or 



{2-nY 

/dz y~ 
— H q (ii, k = — ) D q _> r (z, fi F ; Q) . 



Tq— >7* {k, Q) 



Or 



(5) 



where H q (/i,k = Q + / 1 z) = Tr[S g (//, k = Q + /z){ r y (Q + /z)}} represents the leading power 
short- distance production of a quark of momentum fc M = (Q + /z)n^ and 



D q ^ r (z, fi F ; Q) = j Tr 



(2tt)< 



Q_ 



Tq-ty* (k, Q) 



(6) 



is the quark-to-virtual-photon fragmentation function in terms of its cut-vertex definition 
TTj . As shown in Fig. ||(b), the z 2 5(z — %^ )} in Eq. (||) is the corresponding cut- vertex. 



The cut- vertex for an antiquark-to- virtual-photon fragmentation function is the same as that 
for a quark-to- virtual-photon fragmentation function. 

Similarly, by considering a generic gluon to virtual photon fragmentation process, we 
derive the cut-vertex definition for the gluon-to- virtual-photon fragmentation function with 
the cut- vertex {| d' 11 ' z 2 5(z — as shown in Fig. [l|(c). The tensor d^ v is defined as 



(7) 



From the cut- vertex definitions in Figs. |I|(b) and |T](c) , we derive corresponding opera- 
tor definitions for virtual photon fragmentation functions as follows. By representing the 
diagram in Fig. |I|(b) in terms of quark fields, we have 



D q ^ r (z,/j, F ;Q) = J 



d A k 



(2 



7T 



2 r/ Q + 



4£r 



{2-n)H\k-Q-Y, k x)\{ 



d 3 k 



x 



x 



X 



(2vr) 3 2£ x 



N 



x T7"E Tr 



7 + (0|^(0)| 7 *(Q)X)(X 7 *(Q)|^(0)|0) 



(8) 



dy 

2tt 



-i(Q+/z)y~ 



N 



N 



£Tr 7 + (0|^(0)|7*(g))(7*(Q)l^(^)|0) 



where (l/N)J2iLi with N = 3 indicates the average over the quark color. Similarly, we 
derive the operator definition for the gluon-to-virtual-photon fragmentation function, 
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i N 2 -l 



< 1^7 E (0|^(0)|7*(Q))(7*(Q)l^r(^)|0) • (9) 

a=l 



Both operator definitions in Eqs. (|j) and (H) are in the light-cone gauge. Proper insertion 
of a line integral of the color and electromagnetic potential is needed to make them both 



color and electromagnetic gauge invariant [17,20 



From either the cut- vertex definitions or the operator definitions, we derive the evolution 
equations (or renormalization group equations) for virtual photon fragmentation functions 

+ (?-) E t ^T p c^d(^, a B ) D d ^ r (z\ (i F ; Q) , (10) 
\ZirJ i —f Jz z z 

where c,d — q,q, g. The fragmentation scale up is defined to be the scale where the operators 
of the fragmentation functions are renormalized. We choose the \j? F to be the invariant mass 
of the fragmenting parton at the renormalization point. 

The evolution kernels 7 C ^ 7 * and P c ^d in Eq. ([10|) have the following perturbative expan- 
sions, 

. (z, /iF, a s ;Q) = J2 tSt fa Vf] Q) (^j , (H) 



7c^7 

n=0 

,(n) ^ ( , 

n=0 



W*,*) = E , (12) 



where the renormalization scale dependence is suppressed. In our derivation of Eq. (|H]), we 
assume a strong ordering in partons' invariant mass in the fragmentation sequence, except 
the last step (quark-to- virtual photon). As a result, the evolution kernels P c ^d for the 
homogeneous terms in Eq. flip] ) are the same as the kernels for the DGLAP equation |22 



Although the evolution equations in Eq. (|H]) have the same functional form as that for 
real photon fragmentation functions , the evolution kernels 7 C ^ 7 * are different due to the 
nonvanish mass of the virtual photon. 

In order to calculate the evolution kernels 7 C ^ 7 *, we need to specify the polarization 
vector e%(Q) for the virtual photon state \ j*(Q))- For a unpolarized virtual photon, we need 
only the following polarization tensor, 

p^(Q)= E e T(Q) e x(Q) , (is) 

X=T,L 

where T and L represent the virtual photon's transverse and longitudinal polarization, 
respectively. Although the evolution kernels and the fragmentation functions defined in 
Eqs. (|^) and (|9|) are gauge invariant, the functional form of the polarization tensor P^ U (Q) 
as well as the number of Feynman diagrams contributing to the evolution kernels are gauge 
dependent. 
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In the light-cone gauge, we have the polarization tensor 



Q • n 



(14) 



and have only one Feynman diagram, as shown in Fig. 0(d), which contributes to the lowest 
order quark-to-virtual-photon fragmentation function. We obtain 



Dfl 7 *(z,ii F ;Q) = e\ 



ex. 



2tt 



1 + (1 



In 



Q 2 



(15) 



where e q is the fractional charge for the quark of flavor q, and the fragmentation scale 
A*f > Q 2 / z due to the kinematics. 

For an arbitrary choice of gauge, we need a total of four Feynman diagrams for the lowest 



order quark-to- virtual-photon fragmentation, as shown in Fig. | 



] . The diagrams in Fig. 



contain double "eikonal" lines, and the Feynman rules for the double "eikonal" lines can 
be found in Refs. |T7| , |23|| . These four diagrams in Fig. |2] form a gauge invariant set at this 
order 23,241. Contributions of the three diagrams from Figs. ^](b) to |](d) vanish when they 
are contracted by the light-cone gauge polarization tensor in Eq. fll4"|) . With the covariant 
polarization tensor 



Q 2 



(16) 



or simply P^ U (Q) = —g^, we obtain the same lowest order quark-to- virtual-photon frag- 
mentation function in Eq. (|i~5|) . Since gluon does not directly couple to a photon, we have 



(17) 



By applying fi 2 F d/d[i% to the lowest order virtual photon fragmentation functions, we 
obtain the lowest order evolution kernels 



1 + (1 



0{n 



Q 1 



) 



(0) 



z,fi F ;Q) = 0. 



(18) 



With the lowest order evolution kernels 7^ 7 * in Eq. (|18D and P^i d from Ref. we 
solve the evolution equations in Eq. ( 110 ) for the unpolarized virtual photon fragmentation 
functions, and present the numerical results in Sec. [IV]. 



(0) 



III. POLARIZED VIRTUAL PHOTON FRAGMENTATION FUNCTIONS 

By measuring the momenta of both leptons in Drell-Yan massive lepton-pair production, 
we can determine both invariant mass and polarization state of the virtual photon, which 
decays into the lepton-pair. Therefore, it is meaningful to define the fragmentation functions 
to a polarized virtual photon. In this section, we calculate the evolution kernels for the 
evolution equations of polarized virtual photon fragmentation functions. 
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Since we did not use any constraints due to the virtual photon's polarization when 
we derived the cut-vertex and operator definitions for the virtual photon fragmentation 
functions in last section, the same definitions should be valid for fragmentation functions to 
a polarized virtual photon. Therefore, with a proper choice of polarization vectors €\{Q) for 
a virtual photon of polarization A, we can calculate the fragmentation functions to a virtual 
photon of a specific polarization state. 

In order to specify a polarized virtual photon state, we need to define the photon's 
polarization vector e^(Q). If the photon is real (Q 2 = 0), it has only two transverse po- 
larization states and its longitudinal polarization can be gauged away by an extra gauge 
transformation. When the photon is virtual and its quantum numbers are completely fixed 
by the physical observables, the extra gauge degree of freedom used to remove the longitu- 
dinal polarization at Q 2 = is lost, and therefore, the virtual photon has three polarization 
states. 

With our choice of the frame in which the virtual photon is moving very fast along the 
z-axis, we define the photon's polarization vectors as the following, 

4=i(Q) = (0,1,0,0), 
4 =2 (Q) = (0,0,1,0), 

e£(Q) = -1([Q+ - Q-], 0, 0, [Q+ + Q~\) , (19) 



which are effectively the same as those in the "S-helicity" frame defined in Ref. ||25|| . We 
obtain corresponding polarization tensors, 

PPiQ) = \ E ^t(Q) %(Q) = \ df (20) 
z T=12 z 

with defined in Eq. ([?D and 

Pr(Q)=eT(Q)el(Q) (21) 

for transversely and longitudinally polarized virtual photons, respectively. By summing over 
all polarization states, we should recover the polarization tensor for a unpolarized virtual 
photon, 

P^(Q) = 2 Pf(Q) + Pj?(Q) , (22) 

where the factor of two represents the virtual photon's two transverse polarization states. 

Applying the transverse polarization tensor P^ v to the lowest order Feynman diagrams 
in Fig. 0, we obtain the lowest order fragmentation function for a quark to a transversely 
polarized virtual photon as 



n (0) , 2 {Otem\ 1 ( 1 + (1 ~ Z 



i2 N 




zji 2 Fi 



(23) 



with ftp > Q 2 /z. Again, the lowest order fragmentation function for a gluon to a photon 
vanishes, * (z,jj, F ; Q) = 0. Corresponding evolution kernels are given by 
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(0) 



[z,/i f ;Q) = e, 



1 /l + (l-z) 2 



7, 



(o) 



q 2 
0. 



01 



o 2 



(24) 



,(°) 



7 » in Eq. 



As expected, when Q — > 0, 27 g 
for the real photon fragmentation functions 

Similarly, by applying the longitudinal polarization tensor P£ v in Eq. 
diagrams in Fig. 0, we derive 



4]) reduces to the lowest order evolution kernels 

1|) to the Feynman 



a. 



2tt 



z 



1 



Zfl 2 F/ 



(25) 



with jj, 2 F > Q 2 /z, and ,(z, /j, F ',Q) = 0. Our lowest order fragmentation function to a 
longitudinally polarized virtual photon in Eq. (|25|) is consistent with the result derived by 
Braaten and Lee [p^j . The (1 — z) factor in Eq. (|25| ) is a consequence of the vector interaction 
between the quark and photon. As a consistency check, our lowest order polarized virtual 
photon fragmentation functions in Eqs. (E3f) and fl25|) satisfy 



2 D®. (z, fi F ;Q) + D^> (z, fi F ;Q) = D^> (z, ^ Q) 



j(0) 



(0) 



(26) 



where -D^ 7 * is given in Eq. (|T5|) . From Eq. (p5|), we derive the evolution kernels for longi- 
tudinally polarized virtual photon fragmentation functions 
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(0) 



7 (0) 



[z,v f ;Q) = e 2 q 
'z,ti F ;Q) = 0. 



2 



(27) 



Again, we have 2 7g°l 7 * + 7g-l 7 £ 



7 (0) . 



Since the polarized evolution kernels and the polarized virtual photon fragmentation 
functions are gauge invariant, we can also derive them from the single diagram in Fig. |](d) 
in the light-cone gauge. Substituting = Q^W 1 + Q~n^ into the polarization tensor in the 
light-cone gauge in Eq. (|T4|), we obtain 



Q 2 



(Q 



(28) 



Since the vector used to fix the light-cone gauge does not have transverse component, 
the dT in Eq. (|28|) should still be identified as 2P^(Q). From Eq. @, we obtain the 
polarization tensor for a longitudinally polarized virtual photon 



pf(Q) 



(Q 



(29) 



in the light-cone gauge. Applying this P L M {Q) to the Feynman diagram in Fig. 0(d), we can 
derive the lowest order fragmentation functions to a longitudinally polarized virtual photon. 
As expected, the derived fragmentation functions are the same as those given in Eq. (|25|). 
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Since the fragmentation functions to a polarized and a unpolarized virtual photon share 
the same form of the operator definitions, the evolution equations for polarized virtual pho- 
ton fragmentation functions have the same functional form as that in Eq. (|10"1), except the 
evolution kernels 7 C ^ 7 * are replaced by those in Eq. ([24]) and Eq. ( p7|) for transversely po- 
larized and longitudinally polarized virtual photon, respectively. Since the evolution kernels 
P c ^d are independent of the polarization of the produced virtual photon, they should remain 
the same. 

IV. NUMERICAL RESULTS AND CONCLUSIONS 

In last two sections, we derived the analytical expressions for the lowest order parton- 
to-virtual-photon fragmentation functions with the virtual photon unpolarized as well as 
polarized. The lowest order gluon-to-virtual-photon fragmentation functions vanish while 
the lowest order quark-to- virtual-photon fragmentation functions are given in Eqs. flI5p, 
and (^) for unpolarized, transversely polarized, and longitudinally polarized, respectively. 
Strong interactions in the top blob of the cut-vertex diagrams in Figs. |T](b) and |I](c) can 
change the z as well as the ^f dependence of the lowest order virtual photon fragmentation 
functions. Solving the evolution equations in Eq. is effectively to resum all leading 
logarithmic contributions from the strong interactions. The power corrections from the 
strong interactions 0(Aq CD / fi 2 F ) are not included in the evolution equations. 

In order to solve the evolution equations in Eq. (JT0|) , we need to specify a boundary con- 
dition. Unlike the real photon fragmentation functions, we do not need any nonperturbative 
input distributions if the invariant mass Q ^> Aqcd- From the pure kinematics, we have the 
following boundary condition for solving the evolution equations in Eq. (|10|), 

D f ^ r (z,fi 2 F < Q 2 /z-Q) = (30) 

for all parton flavor / = q,q,g and any polarization of the virtual photon. 

Since the boundary conditions given in Eq. (j3~0l) are the same for all flavors of mass- 
less partons, and the evolution kernels 7 C ^ 7 * in Eq. fllCf ) are the same for a quark q and 
corresponding antiquark q, we have 

D q _> r (z, /ip; Q) = Dq_ T (z, /i F ] Q) (31) 

for all quark flavor q. By neglecting the quark mass difference, the only flavor dependence 
of the evolution kernels 7 C ^ 7 * in Eq. ( |T(]| ) is from quark's fractional charge e q . Therefore, 
quark-to- virtual-photon fragmentation functions are the same for all quark flavors with the 
same fractional charge. 

With the boundary condition given in Eq. (|30|), we can solve the evolution equations 
in Eq. ( |10"D in the moment-space analytically, and then, perform the Mellin transformation 
from the moment-space back to the z-space Q. However, we find that it is easier to solve 
the evolution equations numerically in the z-space directly. 

In Figs. [5](a) and |3](b), we plot the derived lowest order quark-to- virtual-photon frag- 
mentation functions as a function of the momentum fraction z at a fragmentation scale 
(jlf = 10 GeV and \ip = 50 GeV, respectively. We choose e q = 2/3 for the quark's fractional 
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charge, and the virtual photon's invariant mass to be Q = 5 GeV. The unpolarized quark-to- 
virtual-photon fragmentation functions given in Eq. fllBD are represented by the solid lines. 
The transversely and longitudinally polarized virtual photon fragmentation functions given 
in Eqs. fl23|) and ( p5|) are represented by the dashed and dotted lines, respectively. The solid 
lines are equal to twice of the dashed lines plus the dotted lines, which is a consequence of 
Eq. ©. 

From Fig. |3|, we find that longitudinally polarized virtual photon fragmentation functions 
are much larger than the transversely polarized virtual photon fragmentation functions when 
z is small. The small z region corresponds to the region where \j? F is close to the threshold 
Q 2 /z. Near the threshold, the evolution kernels for the transversely polarized fragmentation 
functions 7 g ^ 7 *(0) in Eq. ( f24|) vanish while the kernels for the longitudinally polarized 

fragmentation functions 7^ 7 * in Eq. (p7|) are finite and large. The 7q°i 7 * are actually 
proportional to \ jz when /i F — > Q 2 /z. Therefore, the longitudinally polarized virtual photon 
fragmentation functions dominate the small z or the threshold region. 

On the other hand, the fragmentation functions for a transversely polarized virtual pho- 
ton evolve much faster than the longitudinally polarized fragmentation functions in the large 
z region, as shown in Fig. |3|. This is because the evolution kernels 7^ 7 * for the longitudi- 
nally polarized virtual photon fragmentation functions are power suppressed (proportional 
to 1/fJLp) and also vanish as z — > 1. 

In order to see the effect of resumming the large logarithms from quark and gluon radia- 
tions, we numerically solve the evolution equations in Eq. ( |T0"D with the lowest order evolution 
kernels 7^ 7 . and pJl^j. We plot the comparison between the lowest order (dashed) and the 
QCD evolved (solid) quark-to-virtual-photon fragmentation functions at Q = 5 GeV and 
Hp = 50 GeV in Fig. f|(a). The symbols: U, T, and L represent the fragmentation func- 
tions to a unpolarized, transversely polarized, and longitudinally polarized virtual photon, 
respectively. Like the evolution of parton distributions, QCD evolution tries to enlarge the 
fragmentation functions in the small z region, and to suppress the fragmentation functions 
in the large z region. The difference between the QCD evolved and lowest order fragmen- 
tation functions at Q = 5 GeV is small because of the boundary conditions in Eq. (pOP- 
However, when Q is smaller or \ip is larger, we expect QCD evolution to be much more 
important because of a larger logarithm ln(fi F /Q 2 ). For example, cross sections for virtual 
photon production ware measured by UAl Collaboration at CERN for the virtual photon 
mass Q G [2m^, 2.5] GeV Instead of Q = 5 GeV in Fig. £|(a), we plot the same quark-to- 
virtual-photon fragmentation functions at Q — 1.5 GeV in Fig. |(b). Clearly, QCD evolved 
virtual photon fragmentation functions in Fig. |](b) are enhanced in comparison with those 
in Fig. |](a), particularly, in the small z region. 

Although we do not have the lowest order gluon-to-virtual-photon fragmentation func- 
tions due to 7g_l 7 . = 0, QCD evolution in Eq. ([[(]) can generate the gluon-to-virtual-photon 
fragmentation functions. In Figs. |5](a) and |5|(b), we plot gluon-to-virtual-photon fragmen- 
tation functions at \ip = 10 GeV and fip = 50 GeV, respectively. The virtual photon 
mass is again chosen to be Q = 5 GeV. As shown in Fig. [|, QCD evolution generated 
gluon-to-virtual-photon fragmentation functions at Q = 5 GeV grow very fast when the 
fragmentation scale [ip increases. They are about two orders of magnitude smaller than the 
quark-to- virtual-photon fragmentation functions at fip — 10 GeV, and only one order of 
magnitude smaller at fip = 50 GeV. Therefore, at Q = 1.5 GeV or at a larger value of fip, 
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QCD generated gluon-to-virtual-photon fragmentation functions become more important. 

In conclusion, we have argued that virtual photon fragmentation functions are well- 
defined and physically meaningful. We derive the evolution equations for virtual photon 
fragmentation functions, and show that these fragmentation functions are perturbatively 
calculable. We demonstrate that QCD resummation of the large logarithms caused by quark 
and gluon radiation provides a very important contribution to the fragmentation functions 
when the fragmentation scale fip is large and/or the invariant mass Q is relatively small. 
Just like pion fragmentation functions ||, the virtual photon fragmentation functions derived 
here are universal and can be applied to any processes with massive lepton-pair production. 
In the rest of this paper, we discuss some potential applications of our results. 

Contributions of fragmentation functions D(z) to the physical observables generally de- 
pend on the production of the parent partons at the same z. Since the cross sections for 
producing the parent partons strongly depend on the momenta of the produced partons, 
the role of the fragmentation functions may be different for different regions of the z val- 
ues. For example, for the Drell-Yan production, if the cross section is dominated by the 
small z region, the produced virtual photons will likely be longitudinally polarized. On 
the other hand, if the cross section is dominated by the large z region, the virtual photon 
will be transversely polarized. For a fixed collision energy y/S, the Drell-Yan cross sections 



depend on the z- value from z min = y (Q 2 + Q T )/S [e y + e y ] to 1 with the rapidity y. For 

Q = 5 GeV and y = at \/S = 2 TeV (the new Tevatron energy), z m i n ~ 0.01 and 0.05 
for Qt = 10 and 50 GeV, respectively. It is clear from Fig. |3|(b) that the produced virtual 
photons at Qt ~ 50 GeV are more likely to be transversely polarized. At Qt = 10 GeV, 
the longitudinally polarized fragmentation functions in Fig. |3](a) are much larger than the 
transverse fragmentation functions in small z region, and the 1/z 2 factor in Eq. (||) for the 
convolution over z is also favor for producing a longitudinally polarized virtual photon. But, 
the cross section for producing a parent parton of momentum \zt at small z ~ Qt/^t is a 
very steep falling function of kT, and consequently, it strongly reduces the rate for producing 
a longitudinally polarized virtual photon. In conclusion, when Qt is large, the virtual pho- 
ton produced in Drell-Yan massive lepton-pair production is more likely to be transversely 
polarized. 

Recently, it was found that the J/ip mesons produced at Fermilab Tevatron become more 
longitudinally polarized as the transverse momentum increases |16|| . On the other hand, var- 
ious theoretical calculations predict the J/ip mesons to be more transversely polarized fl27j . 
The virtual photon production (extracted from Drell-Yan massive lepton-pair production) 
at large Qt and small Q 2 has a lot in common with the J /if) production at high Qt- They 
both have two large physical scales: Qt and Q 2 , which is equal to M 2 ,^ in the case of J /if) 

production; and Q\ is much larger than Q 2 . If the collision energy v^S* is large enough 
and the logarithm ln(Q T /Q 2 ) is so important that the resummed fragmentation contribu- 
tions dominate the production cross sections, the virtual photon and J/ip production will 
share the same partonic subprocesses, labeled by the R in Fig. With their respected 
fragmentation functions, both the virtual photon and J /if) production at high Qt are per- 
turbatively calculable PEfl . As discussed above, the virtual photon fragmentation functions 
are completely perturbative, while the parton to J /if) fragmentation functions involve final- 
state nonperturbative interactions. Therefore, only difference between the virtual photon 
and J /ip production at high Qt and large v^S" is the final-state strong interactions during 
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the formation of Z/ip meson. We propose to measure the virtual photon polarization in 
Drell-Yan massive lepton-pair production at large Qt and low Q 2 . The measurements of the 
virtual photon's polarization at high Qt provide not only a good test of QCD perturbation 
theory, but also a reference process to test the models of J/ifi formation. 

As shown in Fig. the fragmentation functions from a parton d to a physical J/i/j can 
be approximated by the fragmentation functions to a virtual gluon of invariant mass Q, 
which immediately decay into a cc-pair, convoluted with a transition from the cc-pair to a 
physical J/tp The fragmentation functions for a parton to a virtual gluon should be 
perturbatively calculable, and their dependence on the polarization should be very similar to 
that of the virtual photon fragmentation functions. But, the produced charm and anticharm 
quark pair of invariant mass Q can in principle radiate gluons and have soft interactions 
with other partons in the collisions. However, due to the heavy quark mass, such final-state 
interactions during the formation or transition from the produced cc-pair to a physical J/ip 
meson are not expected to significantly change the polarization |[27|| . If the formation from 
the cc-pair of invariant mass Q to a physical J/ip meson does not change the polarization, 
one can expect the polarization of the J/i/j mesons produced at high Qt to be similar to the 
polarization of the virtual photon in Drell-Yan massive lepton-pair production at the same 
kinematics. For example, the non-relativistic QCD (NRQCD) model of J/i/j production 
predicts the J/ip mesons to be transversely polarized at large transverse momentum ||27|| , 
which is not consistent with recent Fermilab data ||16|| . Since the leptons do not interact 
strongly once produced, the measurements of the virtual photon polarization in Drell-Yan 
massive lepton-pair production at large transverse momentum should help us to narrow the 
questions about J/ip production. 
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FIGURES 




(c) (d) 



FIG. 1. (a) A generic diagram for a scattering process in which a quark of mo- 
mentum k fragments into a virtual photon of invariant mass Q; (b) the cut- vertex di- 
agram for quark-to-virtual-photon fragmentation function; (c) the cut-vertex diagram for 
gluon-to- virtual-photon fragmentation function; (d) the lowest order cut- vertex diagram for a quark 
to a virtual photon. 
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FIG. 3. The lowest order virtual photon fragmentation functions as a function of z at 
Q = 5 GeV and /Up = 10 GeV (a), and //p = 50 GeV (b). The solid, dashed, and dotted lines are 
for unpolarized, transversely polarized, and longitudinally polarized virtual photons, respectively. 
The inset in Fig. ||(b) shows the z < 0.025 region. 
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FIG. 4. Comparison between the lowest order (dashed) and the QCD evolved (solid) virtual 
photon fragmentation functions as a function of z at fip = 50 GeV and Q = 5 GeV (a) and 
Q = 1.5 GeV (b). 



18 




FIG. 5. QCD generated gluon-to-virtual-photon fragmentation functions as a function of z at 
Q = 5 GeV and \if = 10 GeV (a), and ^f = 50 GeV (b). The solid, dashed, and dotted lines are 
for unpolarized, transversely polarized, and longitudinally polarized virtual photons, respectively. 




FIG. 6. Sketch for Drell-Yan massive lepton-pair and J/tp production via parton fragmentation. 
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